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We investigate the role of entanglement in quantum phase transitions, and show that the success 
of the density matrix renormalization group (DMRG) in understanding such phase transitions is 
due to the way it preserves entanglement under renormalization. We provide a remterpretation of 
the DMRG in terms of the language and tools of quantum information science which allows us to 
rederive the DMRG in a physically transparent way. Motivated by our reinterpretation we suggest 
a modification of the DMRG which manifestly takes account of the entanglement in a quantum 
system. This modified renormalization scheme is shown, in certain special cases, to preserve more 
entanglement in a quantum system than traditional numerical renormalization methods. 
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A fundamental difference between quantum and classi- 
cal physics is the possible existence of nonclassical corre- 
lations between distinct quantum systems. The physical 
property responsible for the nonclassical correlations is 
called entanglement [[jj ||, ||. Recently Q || ||, it has 
been conjectured that entanglement plays a critical role 
in the quantum phase transitions which occur in inter- 
acting quantum lattice systems at zero temperature. 

Unlike ordinary phase transitions, quantum phase 
transitions take place at zero temperature as one or more 
physical parameters of the system are varied. At the crit- 
ical point there are long-range correlations in the system, 
as in a conventional phase transition. However, because 
the system is at zero temperature, and assuming there is 
no ground-state degeneracy, the system must be in a pure 
state. It follows that the correlations which are the prin- 
cipal experimental signature of a quantum phase transi- 
tion are due to long-range entanglement in the ground 
state of the system. 

A variety of methods have been developed to un- 
derstand quantum phase transitions. One of the most 
successful is the density matrix renormalization group 
(DMRG) |l, a numerical technique for finding accu- 
rate approximations to the ground and low-lying ex- 
cited states of interacting quantum lattice systems. The 
DMRG was developed in response to the poor results 
obtained from Wilson's numerical renormalization group 
(RG) |^ when applied to general quantum lattice sys- 
tems. White [|l0| identified the mathematical origin of 
the breakdown of Wilson's RG as the failure of the 
method to properly take into account the effect of bound- 
ary conditions during iterations. 

The focus of this paper is not on proving technical 
results. Rather, the purpose is to provide a novel physi- 
cal picture of the DMRG and quantum phase transitions 
based on the theory of entanglement developed by quan- 
tum information science 0, This physical picture 
provides a fruitful way of interpreting the properties of 
condensed-matter systems, substantially simplifies and 



physically motivates existing derivations of properties of 
the DMRG, and motivates us to propose renormalization 
schemes extending the DMRG. We will explain the physi- 
cal origin of the breakdown of Wilson's RG as a failure to 
take into account the effects of quantum entanglement in 
the ground state of the system, and show that the DMRG 
can be viewed as a natural extension of Wilson's RG to 
take into account the contribution entanglement makes 
to correlations in the ground state. Unfortunately, ow- 
ing to the limited results available on entanglement mea- 
sures, our proposed modifications to the DMRG cannot 
be derived in closed form. For this reason we are forced 
to implement the modified renormalization scheme nu- 
merically for idealised example systems with only small 
numbers of spins. 

The structure of the paper is as follows. First, we 
explain the critical role entanglement plays in any long- 
range correlations present in a quantum system at zero 
temperature. Second, the physical origin of the failure of 
Wilson's RG is identified. The relationship between the 
DMRG and its treatment of entanglement is then elu- 
cidated. Finally, a numerical renormalization scheme is 
presented which has the property of preserving more en- 
tanglement in a quantum system than the DMRG. The 
differences in the entanglement preserved by the vari- 
ous renormalization schemes is illustrated using the ID 
Heisenberg model. 

We now consider the physical origin of the correlations 
which occur in systems exhibiting a quantum phase tran- 
sition. For concreteness, we restrict attention to a lat- 
tice of spin-i particles. If the ground state were a prod- 
uct state, then a simple calculation shows that the spin- 
spin correlation function (erf er^) — (erf) (cr?) is identically 
zero. Thus, if the correlation function is non-zero then 
the ground state must be entangled, and large values of 
the correlation function imply a highly entangled ground 
state. It is well known that for general quantum lattice 
systems the correlation function decays exponentially as 
a function of the separation \i — j\ when the system is far 
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from criticality W]. When the system is at a critical point, 
the correlations decay only as a polynomial function of 
the separation. At this point a fundamental change in 
the ground state has occurred. The physical property 
responsible for these quantum correlations is entangle- 
ment, so that for systems approaching a critical point the 
structure of the entanglement in the ground state must 
undergo a transition. Evidence for such a transition has 
been obtained for the transverse Ising model Q . 

The standard technique used to derive the polynomial 
decay of correlations at the critical point is the renormal- 
ization group Renormalization group techniques as 
applied to quantum field theory and quantum lattice sys- 
tems are now well developed (see, for example, |l3|, [Til , 
for reviews). The renormalization group works by succes- 
sively discarding degrees of freedom in a system until the 
dominant terms may be identified. Each renormalization 
scheme makes a decision about what degrees of freedom 
in a system are important and this decision, in turn, dic- 
tates which class of systems the method is best suited to. 
Prior to 1992 the success of the numerical implementa- 
tion of the renormalization group for quantum systems 
was limited to the Kondo problem |J, and attempts to 
apply Wilson's RG to other quantum lattice systems suf- 
fered from poor convergence and low accuracy. 

Briefly, Wilson's numerical RG procedure works by iso- 
lating a portion A of a quantum lattice system containing 
L sites and exactly diagonalising the Hamiltonian Hl of 
the system restricted to these sites. The subspace W of 
the m lowest eigenvalues and eigenvectors is obtained. 
(Note that m is a parameter we choose. It quantifies 
how many degrees of freedom we keep when renormal- 
izing.) For example, for system A we might pick two 
adjacent spins on the lattice, and then determine the 
two-dimensional subspace W corresponding to the two 
lowest eigenvalues of Hl ■ Once this is done all operators 
are "renormalized" (effectively, averaging out degrees of 
freedom thought to be unimportant), by projecting onto 
the m-dimensional subspace W. The allows us to renor- 
malize the entire lattice, replacing the two spins originally 
in A with a single renormalized spin. On the new renor- 
malized lattice a site is added to A and the interactions 
with the new site are taken into account by constructing 
a new Hamiltonian Hl+i- The exact diagonalisation is 
performed again, keeping the m lowest eigenvalues and 
eigenvectors and the subspace W is updated. The expec- 
tation value of any operator of interest is calculated with 
respect to the ground state of A (which is the approxima- 
tion to the overall ground state), and as the procedure 
is repeated the change in these expectation values are 
noted. Convergence is deemed to have been obtained 
when this difference becomes smaller than some prespec- 
ified value. 

Wilson's RG procedure was expected to work because 
it was assumed that the only degrees of freedom con- 
tributing significantly to the overall ground state would 



be the ground state and low-energy excitations of the lat- 
tice portion A. This assumption turns out not to be the 
case, and Wilson's RG works poorly for most systems. 
White gave a mathematical explanation of why this 
is the case, in terms of the failure of Wilson's RG to prop- 
erly take into account boundary conditions on the system 
being renormalized. Using the tools of quantum informa- 
tion science, we can provide another way to interpret the 
failure of Wilson's RG. The reason the assumption fails 
is that the ground state of an interacting lattice is typi- 
cally entangled, and thus a subsystem of the total lattice 
cannot be assigned a definite state. Indeed, Wilson's RG 
ignores the terms coupling A to the remainder of the 
lattice, which are the very terms responsible for the en- 
tanglement in the ground state causing the correlations 
which signify a quantum phase transition. 

The DMRG improves upon Wilson's scheme by taking 
into account the entanglement between system A and 
the remainder of the lattice. It does this by introducing 
a second subsystem of the lattice, B (the environment), 
in addition to the portion A (now referred to as the prin- 
cipal system) isolated by Wilson's RG procedure. For 
example, if the principal system A is two adjacent spins 
in a one-dimensional lattice, then B might consist of two 
spins, one on either side of A. A and B together are 
referred to as the superblock. Using the physical picture 
based upon entanglement, the DMRG may be described 
as follows. First, construct the superblock Hamiltonian 
and find its ground state. Of course, this ground state 
does not properly describe the entanglement between the 
superblock and the remainder of the lattice. However, 
provided B is chosen appropriately, most of the entan- 
glement between A and the remainder of the lattice will 
in fact be entanglement between A and B. Physically, 
we expect this to happen if B contains all those sites to 
which A most strongly couples. The second step is to 
renormalize A by projecting onto those degrees of free- 
dom of A which contribute most to the entangled ground 
state of the superblock, in a sense made more precise be- 
low. That is, we choose an TO-dimensional subspace W 
of system A, and renormalize the operators for A in a 
similar manner to that described for Wilson's RG. The 
difference between the two approaches is in how W is 
chosen: in Wilson's RG, W is chosen without taking ac- 
count of ground state entanglement, while in the DMRG 
ground state entanglement is critical. This method is it- 
erated until convergence (using the same criteria as for 
Wilson's RG) is obtained. 

Motivated by this physical picture, we can now explain 
the mathematics of the DMRG. Our description is equiv- 
alent to White's |lC| ], but is presented using techniques 
of quantum information science which simplify the proof, 
and highlight the key physical role entanglement plays in 
the DMRG. The central concept we need is the Schmidt 
decomposition p"T[ , one of the main tools used in the 
study of entangled states. The Schmidt decomposition 
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states that a pure state of a bipartite system with 

components A and B may be written as a single sum over 
positive real co-efficients -JpZ, 
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where \u a ) and \v a ) are orthonormal eigenvectors of the 
reduced density matrices p A and p B , respectively. It is 
easy to show that |\&) is entangled if and only if there is 
more than one term in Eq. (|l|). 

Which degrees of freedom should we keep to properly 
take into account the entanglement between A and B? 
One reasonable choice is provided by the DMRG projec- 
tion P, which is defined to be the m-dimensional projec- 
tion on system A that minimizes the functional 



S{P) = m A B-P\*}AB\ 
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The motivation for introducing this functional is to iso- 
late the degrees of freedom which, if discarded, would 
change the state the least. To determine the projection 
which minimizes this functional we write the state of 
the superblock in Schmidt form, \^} =J2 a 
Simple algebra shows that 
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\\V) AB - P\V) AB \ 2 = 1 - tr (pJ2Pa\u a )(u a \j . 
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Thus minimizing the functional in Eq. (||) is equivalent 
to maximizing tr (P^2 a p a \u a )(u a \). The Ky Fan max- 
imum principle |l5[ | tells us that the maximum in this 
expression is achieved by choosing P to be the projector 
onto the space spanned by the \u a ) with the m largest 
values p a . 

Summarizing, the physical origin for the failure of Wil- 
son's RG and the success of the DMRG lies in the way 
each method deals with entanglement. Typically the 
ground state of a lattice system exhibits quantum corre- 
lations. In order to represent the state of an isolated por- 
tion of lattice sites A it is necessary to take into account 
the correlations between A and the rest of the lattice C. 
Due to the nonlocal nature of the entanglement respon- 
sible for the correlations, it is impossible to represent the 
state of AC using only the local degrees of freedom of 
A. Thus Wilson's RG cannot represent the correlations 
of A with C unless a full exact diagonalisation is per- 
formed. The DMRG approximates the correlations of A 
with C by introducing a subsystem B of C in such a way 
that the entanglement between A and B approximates 
the entanglement between A and C. 

We now outline another technique to determine the 
best degrees of freedom to retain during the rcnormaliza- 
tion step that is easily generalized to mixed states. The 
best degrees of freedom to retain may be thought of as 
identifying which m-dimensional subspace will best pre- 
serve the overlap, or fidelity, of the original state with a 



state that has support only on the smaller subspace. We 
write the renormalized state in the m-dimensional sub- 
space in Schmidt form = Y^p=i y/^\ w>3 } l x/? )- The 
quantity we wish to maximize is the fidelity between the 
renormalized and unrenormalized states, (vP^). This 
problem has been considered previously pi 17 where 



it was shown that the maximum fidelity occurs when the 
Schmidt basis for corresponds to the Schmidt basis 
elements of \^>) with the m largest Schmidt coefficients. 
This derivation is therefore equivalent to the projection 
method because the same subspaces are preserved. The 
subspace-fidelity derivation has the added benefit of be- 
ing physically transparent as well as being easily gener- 
alized to mixed states. Indeed, when the state of the 
superblock is mixed, that is, at finite temperature, the 
projection method may be used to show || that the op- 
timal degrees of freedom to retain in the renormalization 
step are the m eigenvectors of the reduced density matrix 
for A with largest eigenvalues. This result may be veri- 
fied with a subspace-fidelity argument where the expres- 
sion for the fidelity of a pair of mixed density matrices is 
F{pi,p 2 ) = tr (^f/pTp^/pl) . 

The DMRG is an excellent numerical scheme for many 
applications because the state after the projection is usu- 
ally close to the original. However, there is strong evi- 
dence that the DMRG docs not reproduce the algebraic 
decay of correlations with respect to separation at the 
critical point jL9), which is the standard signature of a 
quantum phase transition. We conjecture that the origin 
of this failure is due to the fact that while the DMRG 
does preserve some of the entanglement between A and 
the remainder of the lattice, it does not preserve the max- 
imal amount of entanglement. To make this intuition 
more precise, we need a good quantitative measure of 
how much entanglement is present in a quantum state. 
For pure states such a measure of entanglement between 
A and B is the von Neumann entropy of the individual 
systems A and B alone, S(A) = —trpA logp A = S(B) M. 
We conjecture that the reason the DMRG does not re- 
produce the algebraic decay of correlations at the critical 
point is because it minimizes the difference (^|) rather 
that the difference between the entanglement of AB and 
the entanglement of AB after the projection. 

This conjecture motivates the construction of a projec- 
tion preserving more entanglement between the system 
and environment than the DMRG. For the case where the 
superblock is in a pure state we wish to determine which 
projection P on A when applied to the joint state | v I')ab 
of AB will preserve the maximal amount of entanglement 
between A and B. This is equivalent to maximizing 



T = S(P PA P/tr(Pp A )), 
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over all m-dimensional projections P on A. To do this 
maximization numerically, we introduce a basis {|</>i)} 
for A and write P = \4>i){4>i\ where the sum runs 
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over some subset of the basis. The projection P is deter- 
mined by varying the directions of the basis vectors until 
the quantity in (Q) is maximized. Numerically this is 
done by introducing a Lagrange multiplier for each basis 
vector to enforce orthonormality. We refer to the projec- 
tion P constructed by the procedure outlined here as the 
entanglement maximization projection (EMP). 

For nonzero temperatures the state of the superblock 
Pab is mixed. The derivation of the EMP outlined in 
the previous paragraph requires a measure of the entan- 
glement of A with B for mixed states in order to maxi- 
mize (||). For our calculations we use the best- understood 
measure of mixed-state entanglement, the entanglement 
of formation E(A : B) || of the joint state pab- The 
construction of the EMP is otherwise is unchanged. Un- 
fortunately, no simple formula for the entanglement of 
formation has yet been obtained. For this reason we are 
only able to calculate the EMP for systems with small 
numbers of spins. When a simple formula for a good 
mixed-state entanglement measure is obtained it may be 
possible to derive a general formula for the EMP. 

The EMP has the property that entanglement in the 
system is preserved maximally. The EMP does not min- 
imize the functional Eq. (|J) so that, in general, the state 
is quite different after the projection. For this reason the 
EMP is not appropriate for most calculations. However, 
there is one situation where we expect that the EMP may 
be useful, that is, for calculating properties of a system at 
a quantum critical point. At this point the correlations, 
which are due to entanglement, are strongest. We con- 
jecture that the EMP (unlike the DMRG) may be able 
to reproduce the algebraic decay of correlations at the 
critical point, but have not yet been able to verify this 
analytically or numerically. 

The difference between the entanglement preserved by 
the EMP and the DMRG projection may be seen us- 
ing the example of the ID antiferromagnetic Heisenberg 
model H = J/2J2 (1 ?) ct 1 ■ ^ (See [| || || || 



for other computations of the entanglement in condensed 
matter systems.) As we noted earlier, for computational 
ease we chose the model to be on four sites. The EMP 
was calculated by optimizing over the entanglement of 
formation E(A : B) of the projected system using the 
formula of Wootters p5j. The entanglement of forma- 
tion as a function of temperature is shown in Fig. [I]. It is 
clear from Fig. |l| that the EMP better preserves the en- 
tanglement present in the system than does the DMRG, 
as expected. Similar results have been obtained for the 
Ising model in a transverse field, and we expect that this 
is typical of strongly interacting quantum lattice systems. 

In this paper we have shown that entanglement plays 
a central role in quantum phase transitions. We have 
shown that the success of the DMRG — one of the prin- 
cipal numerical techniques used to solve quantum lattice 
systems near criticality — is due to the way it takes ac- 
count of the entanglement in such systems. We have also 



1.4 
1.2 
1 

0.8 
0.6 
0.4 
0.2 






Lower bound 


\ 

\ 

» — . 

1 \ 
I \ 
\ \ 


EMP 

DMRG 


1 \ 
1 \ 
\ \ 

~ v \ \ 

N <i \ 




\N \ 
'> \ \ 




i s \ 




1 v 









3 

kT 



FIG. 1: Sites 1, 2 in the 4-site Heisenberg model were renor- 
malized to a single spin via the DMRG and EMP projections 
(similarly for sites 3, 4) and the entanglement between the two 
renormalized spins plotted. Also shown is a lower bound on 
the actual entanglement present between A and B obtained 
from the entropy-entanglement inequality (chapter 8 in ]B|, 
and |^|). The lower bound is expected to be a good ap- 
proximation only when kT < 1, where we see that the EMP 
provides a much better approximation to the entanglement 
present in the system than does the DMRG. 



proposed a modification of the DMRG which preserves 
more entanglement during the renormalization step, and 
which we conjecture may be capable of reproducing the 
polynomial decay of correlations near the critical point. 
We believe our results show that the techniques of quan- 
tum information science can be fruitfully applied to ob- 
tain fundamental insights into the properties of many- 
body systems. 
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